Electric Quadrupole Moments of Metastable States of Ca + , Sr + , and Ba 4 



OO 
O 

o 

(N 






Dansha Jiang and Bindiya Arora 
Department of Physics and Astronomy, University of Delaware, Newark, DE 19716-2570, USA 

M. S. SafronovE0 

Physique des Interactions loniques et Moleculaires (CNRS UMR 6633), 

Universite de Provence, Centre de Saint-Jerome, 

Case C21, F13 391 Marseille Cedex 20, France and 

Department of Physics and Astronomy, University of Delaware, Newark, DE 19716-2570, USA^ 

Electric quadrupole moments of the metastable nd 3 / 2 and nd 5 / 2 states of Ca + , Sr + , and Ba + 
are calculated using the relativistic all-order method including all single, double, and partial triple 
excitations of the Dirac-Hartree-Fock wave function to provide recommended values for the cases 
where no experimental data are available. The contributions of all non-linear single and double 
terms are also calculated for the case of Ca + for comparison of our approach with the CCSD(T) 
results. The third-order many body perturbation theory is used to evaluate contributions of high 
partial waves and the Breit interaction. The remaining omitted correlation corrections are estimated 
as well. Extensive study of the uncertainty of our calculations is carried out to establish accuracy 
of our recommended values to be 0.5% - 1% depending on the particular ion. Comprehensive 
comparison of our results with other theoretical values and experiment is carried out. Our result 
for the quadrupole moment of the 3d 5 / 2 state of Ca + ion, 1.849(17) ea\, is in agreement with the 
most precise recent measurement 1.83(1) ea\ by Roos et al. [Nature 443, 316 (2006)]. 
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I. INTRODUCTION 



Frequency standards based on optical transitions of 
trapped ions have the potential to reach a systematic 
fractional uncertainty on the order of 10 -18 [l|. The 
ability to develop more precise optical frequency stan- 
dards will open ways to improve Global Positioning Sys- 
tem (GPS) measurements and tracking of deep-space 
probes, perform more accurate measurements of the 
physical constants and tests of fundamental physics such 
as searches for nonlinearity of quantum mechanics, grav- 
itational waves, etc. Some of the promising candidates 
for such ultra-high-precision frequency standards with 
trapped ions are 27 A1+ @,|jj, 199 Hg+ ||J|, 171 Yb+ 0,0, 
87 Sr+ g||, 43 Ca+ [I|, 115 In+ [ll|, and 137 Ba+. One of 
the largest sources of systematic errors in such frequency 
standards with monovalent ions is due to interaction of 
the quadrupole moments of metastable states with stray 
electric field gradients [lj, [l3| . The electric quadrupole 
moments of the metastable states are hard to calculate 
accurately even for simplest monovalent systems owing 
to large correlation corrections (over 30% for Ca + ). Rel- 
ativistic configuration interaction (RCI) method with a 
multiconfiguration Dirac-Fock orbital basis was used by 
Itano [14| to calculate relevant quadrupole moments in 
Ca+, Sr+, Ba+, Yb+, Hg+, and Au. The RCI results 
agreed with available measurement within 10%. The rel- 
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ativistic coupled-cluster calculations of quadrupole mo- 
ments of metastable nd states were carried out by Sur et 
al. [H for Ca+, Sr+, and Ba+ and by Sahoo [H for 
Ba + . These calculations yielded results 5% and 13% 
higher than the recent measurements of Ca + 17] and 
uadrupole moments, respectively. Mitroy and 
[201 ] calculated the quadrupole moments of the 
3^5/2 state in Ca + and the id 5 1 2 state in Sr + by diagonal- 
izing a semiempirical Hamiltonian in a large dimension 
single electron basis. Their values are in good agree- 
ment with the experiment. However, they have noted 
that their particular definition of the polarization po- 
tential may lead to a possible problem with the accu- 
racy of properties of these nd states calculated using this 
method [19|, [201 . Large differences between theoretical 
calculations and experimental values, especially the 5% 
discrepancy between the recent precise measurement of 
the 3d 5 i 2 state quadrupole moment of Ca + by Roos et 
al. [17| and accurate coupled-cluster single-double partial 
triple [CCSD(T)] value from [15J . and consecutive need 
for comprehensive analysis of the theoretical uncertain- 
ties have in part motivated this work. 

In this paper, we present the relativistic all-order calcu- 
lations of the electric quadrupole moments of the nd 3 / 2 
and nd 5 / 2 states of Ca + , Sr + , and Ba + ions. The rel- 
ativistic all-order method is one of the most accurate 
methods used for the calculation of atomic properties of 
monovalent systems (see Ref. [211 ] for a review and refer- 
ences therein) . The lifetimes of the 3d levels in Ca + cal- 
culated in this approach and estimated to be accurate to 
1% were found to be in agreement with the high-precision 
experiment [22j. The calculation of the nd quadrupole 
moments is very similar to the calculation of the nd life- 



times so similar accuracy is expected. The long lifetimes 
of the metastable nd states of these ions also make these 
systems well suited for the study of quantum information 
processing and quantum simulation [23, [2J] . 

The atomic properties of Ba + are also of particular 
interest owing to the prospects of studying the parity 
nonconservation with a single trapped ion [251 ] . Progress 
on the related spectroscopy with a single Ba + ion is re- 
ported in [26l . \2J\ , and precision measurements of light 
shifts in a single trapped Ba + ion have been reported in 

m. 

Another motivation for this work is an opportunity to 
evaluate the importance of the non-linear terms as well 
as triple and higher excitations in the coupled-cluster ap- 
proach. It has been indicated in Rcfs. 29, 30] that non- 
linear terms may be relatively large and significantly can- 
cel with triple and higher-excitation terms that are not 
included in the CCSD(T) of Ref. p]| or all-order single- 
double partial triple (SDpT) approaches [21[. In this 
work, we include all non-linear terms at single-double 
(SD) level and evaluate triple and higher-excitation cor- 
rections beyond CCSD(T) or SDpT treatments for Ca + . 
Our results demonstrate significant cancelation between 
these terms. We note that it is the first time that such 
calculation has been carried out for any of the atomic 
properties. Our results in part explain the discrepancy 
of previous high-precision calculation with the Ca + ex- 
periment and represent the most complete calculation to 
date. We also present detailed analysis of the uncertainty 
of our calculations. 



II. METHOD 

The electric quadrupole moment 0(7 J) of an atom in 
electronic state I7 J) is defined as the diagonal matrix el- 
ement of the q = component of the electric quadrupole 
operator Q in a spherical basis 



0( 7 J) = (VfrJMj) |Q | ^(jJMj)} 



(1) 



with the magnetic quantum number Mj taken to be 
equal to its maximum value, Mj = J 14]. Applying 
the Wigner-Eckart theorem and using analytical expres- 
sion for the relevant 3-j coefficient [3l[ , allows to express 
the quadrupole moment via the reduced matrix element 
of the quadrupole operator as 



6(7.7) = 



(2.7)1 



y/(2J-2)\(2J + 3)\ 



<tt( 7 J) ||Q|| tt( 7 J)> , (2) 



where electric quadrupole operator Q is represented in 
second quantization as a one-body operator 



Q = ^2<Hj4 a J 



(3) 



'■j 



Here, a\ and a* are the creation and annihilation opera- 
tors. 



In the coupled-cluster method, the exact many-body 
wave function $(7 J) is represented in the form [32j 

\*)=exp(S)\^), (4) 

where "I^ )) is the lowest-order wave function. Wc 
have omitted indexes (7 J) in this equation and formu- 
las below for convenience. The operator S for an N- 
electron atom consists of "cluster" contributions from 
one-electron, two-electron, ■ • ■ , N-electron excitations of 
the lowest-order wave function |\Ij( )): 



S — Si + 02 



Sn- 



(5) 



The expansion of the exponential in Eq. ^ in terms of 
the n-body excitations S n gives 

I*) = |l + Si + S 2 + ^Sf + 5i5 2 + I^| + ... j |tf(°)). 

.. (6) 
In the linearized coupled-cluster method, only linear 

terms are considered, and the wave function takes the 

form 



I*) 



{1 + Si + s 2 + s 3 



S N }\*V>). (7) 



We note that the contributions from the nonlinear terms 
are expected to be relatively small, but the computa- 
tional complexity and time increases sig nificantly with 
their addition in the present approach 33] . The relativis- 
tic all-order single-double (SD) method is the linearized 
coupled-cluster method restricted to single and double 
excitations only, with the wave function given by 



I* 



SD/ 



{1 + Si + S 2 }|^°)) 



1 + ^2 P->naain a a +^2 



PmvQ j rnP"v 



Et f 
PrrwiabQ'jnQ'n^'b&a 



£. 



^mnva ^rn ^n ^a ®"V 



mnab 



I* 



(8) 



(o)\ 



where we take frozen-core Dirac-Hartree-Fock (DHF) 
wave function to be the lowest-order wave function l^ )}. 
The indices m and n designate excited states while in- 
dices a and b designate core states; the index v labels 
the valence electron. The equations for the single excita- 
tion coefficients p ma , Pmv, double excitation coefficients 
Pmnabi Pmnva, and the corresponding correlation core and 
valence energies SE core , 5E V are solved iteratively in a 
finite basis set. The finite basis set used in our calcula- 
tions consists of single-particle orbitals which are linear 
combinations of B-splines [341 ] constrained to a spherical 
cavity. 

The all-order SDpT method is an extension of the SD 
method in which valence part of the linear triple excita- 
tion term S3 is added to the wave function: 

|*sd p t) = |*sd) 

+ J^ Pm n ™a&aL a n a r a & a a a i>l* (0) ), (9) 



where the |\&sd) is given by Eq. ([8]). The dominant part 
of S3 is treated perturbatively, i.e. its effect on the va- 
lence energies SE V and single excitation coefficients p mv is 
calculated, but the equations for the triple excitation co- 
efficients Pmnrvab are not iterated. A detailed des crip tion 
of the SD and SDpT methods is given in Refs. [HHSSl- 

We carry out both SD and SDpT calculations in this 
work to establish the size of the triple corrections in the 
perturbative approach. The coupled-cluster CCSD(T) 
method used in calculation of the quadrupole moments 
in Refs. [la, [lj| also includes the triple excitations per- 
turbatively even though the particular terms that are 
considered may somewhat differ. 

In this work, we also carry out the all-order calculation 
that includes all non-linear terms that arise from single 
and double excitation terms Si and S 2 for the case of 
Ca + . There are only six of such terms that can contribute 
to the equations for single and double excitation coef- 
ficients, and the complete coupled-cluster single-double 
(CCSD) wave function is then written as 

I^ccsd) = exp(£i + S 2 )\¥V) = |* SD ) + (10) 

[^sf + Sl s 2 + \sl + l -sf + l -sfs 2 + !#} |*<°>>, 

where the |^sd) is given by Eq. IjHJ). The complete formu- 
las for the CCSD equations are given in Ref. [33]. Our 
approach allows us to explicitly calculate the contribu- 
tion of the non-linear terms to the quadrupole moments 
as the difference of the results obtained in the CCSD and 
SD approaches. 

The matrix element of any one-body operator Z in the 
all-order method is obtained as 



(Vv\Z\y w ) 



y/{* v \* v ){V w \* w ) 
For non-scalar operators, this expression becomes 

7 -Zval 

^vw — n rr ri 

^/(1 + N v )(l + N w ) 



(11) 



(12) 



where the expression for the numerator of Eq. (|12|) de- 
rived with I^sd) wave function consists of the sum of the 
DHF matrix element z wv and twenty other terms Z^ k > , 
k = a ■ ■ -t. These terms and the normalization terms 
N v are linear or quadratic functions of the excitation co- 
efficients p ma , Pmv, Pmnab, and p m nva- The complete 
expression for the matrix elements can be found in [37| . 
The same expression [Eq. (ll2[) ] for the matrix elements is 
used in all calculations in this work. 

We carry out three different ah initio calculations of 
the quadrupole matrix elements. In the first one, all 
excitation coefficients are obtained in the SD approach 
[Eq. J5J], in the second one they are obtained in the 
SDpT approach [Eq. ©], and in the third one (carried 
out for Ca + ) the excitation coefficients are obtained in 
the CCSD approach [Eq. fTU)) ]. We refer to these results 



as SD, SDpT, and CCSD values in the text and tables 
below. 

While the numerator of the Eq. (|12[) contains twenty 
correlation terms, only one term is overwhelmingly domi- 
nant for the quadrupole moments considered in this work, 
contributing over 90% of the total correlation correction. 
Following the notation of Ref. (37( , this is the term Z^ c ' 
that is equal in the case of the diagonal quadrupole ma- 
trix element to 



Z (c) 



*£■ 



(13) 



where sum over m ranges over all excited basis set 
states. The lowest-order DHF matrix elements q^ of the 
quadrupole operator are given by 

/>oo 

Qij = <*HC (2) ||j) / r 2 ( 9i (r) 9j (r) + h{r)f 3 (r)) dr, 



(14) 
where C' 2 ' is the normalized spherical harmonic of rank 
2 and <?,-, fa are large and small components of the Dirac 
wave function, respectively. The p mv are single valence 
excitation coefficients calculated in either SD [Eq. ©], 
SDpT [Eq. ©], or CCSD [Eq. JTDJ)] approximations as 
described above. Therefore, evaluation of the omitted 
higher-order corrections to p mv provides an estimate of 
the dominant part of the missing contributions in each 
approximation. These excitation coefficients are closely 
related to the correlation energy SE V . If we introduce 
the self-energy operator (also referred to as correlation 
potential in some works) E m „ as 



(e« 



8E V ) p T , 



(15) 



where ti is the DHF energy of the state i, then the corre- 
lation energy would correspond to the diagonal term £ TO . 
Therefore, the omitted correlation correction can be esti- 
mated by adjusting the single-excitation coefficients p mv 
to the experimentally well-known value of the valence 
correlation energy, and then re-calculating the matrix el- 
ements using Eq. (TT2")) with the modified coefficients [35f 



_5E^_ 
'SE v heoiy ' 



(16) 



The 5E° xpt is defined as the experimental energy [38J mi- 
nus the lowest order DHF energy e v . The theoretical cor- 
relation energy is somewhat different in the SD, SDpT, 
and CCSD approaches. Therefore, this scaling procedure 
has to be conducted separately for each of these three cal- 
culations with SEl hcor y taken to be <5^ D , 5£f DpT , and 
6Ey SD , respectively. We refer to the results of these 
calculations as SD SC , SDpT sc , and CCSD SC values. 

Before discussing the final results of our calculations, 
we describe the calculation of two other corrections that 
need to be accounted for in the ab initio SD, SDpT, and 
CCSD calculations. Any sum over the excited states in 



TABLE I: Contributions of high partial waves and Breit interaction to the electric quadrupole moments of Ca + calculated using 
third-order many-body perturbation theory; l m ax is the highest number of partial waves included in the particular calculation. 
All values are given in atomic units. 



State 


i>max — 


''max — c5 


''max — -LU 


''max — J-Z 


Z = 7. . . 12 


Breit 


3rf 3 / 2 
3^5/2 


1.134 
1.628 


1.127 
1.617 


1.124 
1.614 


1.123 
1.612 


-0.011 
-0.015 


-0.001 
-0.003 



either the calculation of the excitation coefficients or ma- 
trix elements using the Eq. (|12p involve the sum over 
the principal quantum number, calculated essentially ex- 
actly, and the sum over the partial waves, that needs to 
be truncated after some value l ma x (see the sum over m 
in Eq. (fT3|) for an example). In all of our all-order calcula- 
tions, we chose l ma ,x = 6. We find that the contributions 
from higher partial waves are small but significant and 
should not be omitted at the present level of accuracy. 
The size of the contribution of the higher partial waves 
may also shed some light on the disagreement of some 
previous calculations with experiment. 

To evaluate this contribution, we first carried out a 
third-order many-body perturbation theory (MBPT) cal- 
culation with the same basis set and l ma x as the all-order 
calculations, and then performed the same calculation 
with larger basis set and larger l max . A detailed de- 
scription of the third-order MBPT method is given in 
Ref. [H. 

The results of the third-order calculation with the in- 
creasing values of l ma x for the quadrupole moments of 
Ca + are given in Table HI While the total contribution 
of the ' = 7, 8 partial waves is rather substantial, 0.6%, 
contributions of even higher partial waves are small. We 
truncated the sum after / = 12, with the expected un- 
certainty of this truncation being well below our final 
accuracy. The difference between the MBPT calculation 
with lmax = 6 and l max = 12 is taken to be the correction 
for the contribution of higher partial waves and is added 
to the ab initio all-order results. 

We have also evaluated the total contribution of the 
I = 5, 6 partial waves to establish its size, and found it 
to be 3%. Moreover, the inclusion of larger number of 
partial waves reduces the values of the quadruple mo- 
ments since term Z^ c ' contributes with the sign opposite 
to that of the lowest-order value. The inclusion of larger 
number of partial waves increases the absolute value of 
the correlation correction leading to lower total final val- 
ues. Therefore, omitting contributions of higher partial 
waves or exclusion of such orbitals from the basis set in 
other calculations may result in an overestimation of the 
quadrupole moments by a few percent. 

We also investigated the effect of the Breit interaction 
which arises due to exchange of a virtual photon between 
atomic electrons and can be written as 

Bu = a.i ■ a, + - — [ a t ■ (Xj - (q, • r i3 ) (<Xj ■ r y )] , 



2n 



(17) 



where ctj are the Dirac matrices. This correction in- 
cludes instantaneous magnetic interaction between Dirac 
currents (the first term) and the retardation correction 
to the electric interaction (the second term). In order 
to calculate the Breit correction to the quadrupole ma- 
trix elements, we modify the generation of B-spline basis 
set to intrinsically include the Breit interaction on the 
same footing as the Coulomb interaction, and repeat the 
third-order calculation with the modified basis set. The 
difference between the new values and the original third- 
order calculation is taken to be the correction due to 
Breit interaction. This contribution is listed in the last 
column of Table [JJ In second quantization, Breit interac- 
tion operator in a normal form is separated into a one- 
body part and a two-body part [4CJ . The two-body Breit 
contribution is omitted in our approach. The total Breit 
corrections are small and are below the estimated uncer- 
tainty of our theoretical values. Therefore, the possible 
uncertainty introduced by the omission of the two-body 
Breit correction is negligible. In fact, we find that most 
of the Breit correction arises at the DHF level. 



III. RESULTS AND DISCUSSION 

The results of our calculations of the quadrupole mo- 
ments of the metastable nd s / 2 and nd 5 / 2 states of Ca + , 
Sr + , and Ba + ions are summarized in Table QTJ where 
we list the lowest order DHF, third-order MBPT, all- 
order SD and SDpT ab initio, and corresponding all-order 
scaled values calculated as described in Section [Til In 
the case of Ca + , we also list the results of our CCSD 
and scaled CCSD calculations. The ab initio values con- 
tain the corrections for the higher partial wave contri- 
butions and Breit interaction. These corrections do not 
need to be included into the scaled results as that will 
lead to double counting of these effects. We take the 
scaled SD numbers as the final values based on the com- 
parisons of similar calculations in alkali-metal atoms with 
experiment (see Refs. [4l], [42|, [43|, [44], J45| and references 
therein). 

We take the maximum difference between the final val- 
ues and the SDpT ab initio, SDpT scaled [SDpT sc ], and 
CCSD scaled [CCSD SC ] values to be the uncertainty of 
the dominant contribution. We assume that any remain- 
ing uncertainty does not exceed the uncertainty of the 
dominant term and take it to be equal to the uncertainty 
in the dominant term evaluated as described above. The 
two uncertainties are added in quadrature to obtain the 



TABLE II: Electric-quadrupole moments of Ca + , Sr + , and Ba + calculated using different approximations: Dirac-Hartree-Fock 
(DHF), third-order many-body perturbation theory (MBPT), single-double all-order method (SD), and single-double all-order 
method including partial triple excitation contributions (SDpT), label "sc" represent the corresponding scaled values. The 
results of the full single-double couple-cluster calculation for Ca + are listed in row labeled "CCSD" ; the corresponding scaled 
values are listed in the row "CCSD 3C ". All values are given in atomic units. 



Ion 


State 


DHF 


MBPT 


SD 


SDpT 


SD SC 


SDpT sc 


CCSD 


CCSD SC 


Final 


Ca+ 


3^3/2 


1.712 


1.122 


1.245 


1.282 


1.289 


1.281 


1.271 


1.292 


1.289(11) 




3(4/2 


2.451 


1.610 


1.785 


1.837 


1.849 


1.836 


1.822 


1.851 


1.849(17) 


Sr+ 


4d 3 / 2 


2.469 


1.876 


1.987 


2.021 


2.029 


2.020 






2.029(12) 




4d 5 /2 


3.559 


2.721 


2.876 


2.922 


2.935 


2.923 






2.935(17) 


Ba+ 


5^3/2 


2.732 


2.086 


2.217 


2.260 


2.256 


2.248 






2.256(11) 




5d 5 /2 


3.994 


3.087 


3.263 


3.323 


3.319 


3.308 






3.319(15) 



TABLE III: Comparison of the present results for electric quadrupole moments in Ca + , Sr 4 
and experiments. All values are in atomic units. 



and Ba with other calculations 



Ion 


State 


Present 


Ref. [14] 


Ref. [15] 


Other 


Expt. 


Ca+ 


3^3/2 


1.289(11) 


1.338 


1.338 








3<4/2 


1.849(17) 


1.917 


1.916 


1.819" 


1.83(1) 6 


Sr+ 


4d 3 / 2 


2.029(12) 


2.107 


2.12 








4d 5 /2 


2.935(17) 


3.048 


2.94(7) 


2.840 c 


2.6(3) d 


Ba+ 


5^3/2 


2.256(11) 


2.297 


2.309 


2.315 e 






5rf 5 /2 


3.319(15) 


3.379 


3.382 


3.382 e 





"Reference 
'Reference 
c Reference 
d Reference 
'Reference 



19] 



CCSD(T) 



final estimate of the uncertainty of our values. 

We make several conclusions from our results (all the 
% numbers are given for Ca + , but the general trends are 
the same for all ions considered in this work) : 

1. The triple contributions included in the perturba- 
tive approach contribute about 3% and increase the 
values of the quadrupole moments. 

2. The non- linear terms contribute about 2% and also 
increase the values of the quadrupole moments. 

3. While the SD, SDpT, and CCSD results vary by 
a few percent, the addition of the estimated omit- 
ted correlation correction carried out according to 
Eq. (fll)|) , brings all these results to very close agree- 
ment providing additional validation of this proce- 
dure. 

4. The linearized SD SC and complete coupled-cluster 
CCSD SC scaled results are nearly exactly the same, 
with the differences being well below our estimated 
uncertainty. Therefore, we found it unnecessary to 
carry out CCSD calculations for Sr + and Ba + . 

5. We confirm that non- linear terms strongly cancel 
with the triple and higher-excitation contributions 
not included in the perturbative approach. As a re- 
sult, CCSD(T) method used in Refs. [U, \M that 



includes both non-linear terms and triple excita- 
tions in the perturbative approach only, is expected 
to yield results a few percent higher than the ex- 
perimental values. 



In Table IIII[ we compare our final values with other 
calculations and available experimental results. We note 
that our calculation is the most complete one at the 
present time. 

The J-independent moments, i.e. the values with the 
all angular factors divided out, can be obtained by multi- 
plying our results in Table ITU by 5 and 7/2 for the nd^/2 
and nd 5 /2 states, respectively, according to Eq. UJ. 

Our values are systematically lower than the results 
of relativistic configuration interaction (RCI) calculation 
carried out with a multiconfiguration Dirac-Fock orbital 
basis by Itano [14j . As we noted above, high partial 
waves ( I > 4) contribute significantly (about 4%) to the 
quadrupole moments and reduce the values. Therefore, 
the restriction of the excitations to mostly low-/ orbitals 
in Ref. [lj| is expected to lead to higher values in RCI 
calculations. The relativistic coupled-cluster CCSD(T) 
results by Sur et al. [15( for Ca + , Sr + , and Ba + and by 
Sahoo [16| for Ba + are also systematically lower than our 
values, with the exception of the 4^/2 Sr + quadrupole 
moment, which is in agreement with our value. It is un- 
clear why this one value compares differently Since this 
particular value was the focus of the work [15| , perhaps it 



was treated differently from the other cases. As we noted 
above, we expect the CCSD(T) results of Refs. [H, OJ] 
to be a few percent too high owing to the cancelation 
of the non-linear terms and higher-excitation terms not 
included in CCSD(T) approach. Another possible issue 
is the treatment of the high partial wave contributions. 
While the tests of various basis sets were conducted in 
Ref. 15], it is not stated how high partial waves were con- 
sidered. We note that the im plem entation of the coupled- 
cluster method in Refs. [la, [l6[ is significantly different 
from ours and is more closely related to the quantum 
chemistry calculations. 

The results of Mitroy and Zhang [HJ, [20| calculated 
by diagonalizing a semiempirical Hamiltonian in a large 
dimension single electron basis are in good agreement 
with experiment. Our analysis of the correlation correc- 
tion is consistent with such results. We demonstrated in 
Section |TT] that the dominant part of the correlation cor- 
rection to quadrupolc matrix elements comes from the 
term containing essentially the correlation potential E„ m 
that is closely related to the correlation energy. Since 
the cutoff function in the semiempirical potential used in 
Ref. [PJ, [20| is adjusted to reproduce experimental bind- 
ing energies, it appears to be a good representation for 
this application. 

Our result for the 3d 5 / 2 Ca + quadrupolc moment 
1.849(17) eag agrees within the quoted uncertainties with 
the recent high- precision measurement 1.83(1) eag re- 
ported by Roos et al. in Ref. [17| that was carried out 
using a decoherence-free subspace with specially designed 
entangled states of trapped ions. We also verify (by vary- 
ing the nuclear parameter) that our value is not depended 



on the particular isotope within our accuracy. Our result 
for the Sr+ 4d 5 / 2 quadrupole moment 2.935(17) eag is 
in good agreement (just outside of the upper Ict bound) 
with experimental value 2.6(3) eag by Barwood et al. [181 ] 
measured with a single laser-cooled ion confined in an end 
cap trap with variable dc quadrupole potential. 



IV. CONCLUSION 

In summary, we performed a relativistic coupled clus- 
ter calculation of the electric quadrupole moments for 
the nd 3 / 2 and nd 5 / 2 states of Ca + , Sr + , and Ba + ions. 
Our analysis of various contributions in part explains the 
discrepancy of previous high-precision theory with exper- 
iment. We also present detailed evaluation of the uncer- 
tainty of our results and provide recommended values for 
the cases where no precision experiments are available. 
Our result for the quadrupole moment of the 3d 5 / 2 state 
of Ca + ion 1.849(17) ea\ is in agreement with the recent 
measurement 1.83(1) eag by Roos et al. [17l ]. 
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